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ABSTRACT 

Context. The physical conditions of the solar photosphere change on very small spatial scales both horizontally and vertically Such 
a complexity may pose a serious obstacle to the accurate determination of solar magnetic fields. 

Aims. We examine the applicability of Milne-Eddington (ME) inversions to high spatial resolution observations of the quiet Sun. Our 
aim is to understand the connection between the ME inferences and the actual stratifications of the atmospheric parameters. 
Methods. We use magnetoconvection simulations of the solar surface to synthesize asymmetric Stokes profiles such as those observed 
in the quiet Sun. We then invert the profiles with the ME approximation. We perform an empirical analysis of the heights of formation 
of ME measurements and analyze the uncertainties brought about by the ME approximation. We also investigate the quality of the fits 
and their relationship with the model stratifications. 

Results. The atmospheric parameters derived from ME inversions of high- spatial resolution profiles are reasonably accurate and can 
be used for statistical analyses of solar magnetic fields, even if the fit is not always good. We also show that the ME inferences 
cannot be assigned to a specific atmospheric layer: diff'erent parameters sample diff'erent ranges of optical depths, and even the same 
parameter may trace diff'erent layers depending on the physical conditions of the atmosphere. Despite this variability, ME inversions 
tend to probe deeper layers in granules as compared with intergranular lanes. 

Key words. Sun: magnetic fields - Sun: photosphere - Instrumentation: high angular resolution 



1. Introduction 

The solar spectrum carries information about the properties of 
our star. In general, a broad range of atmospheric layers con- 
tribute to the shape of the spectral lines, making it difficult to 
extract this information directly. Both the measurement process 
and the method of analysis introduce uncertainties in the phys- 
ical quantities retrieved from the observations. Sources of er- 
ror are photon noise and instrumental effects like limited spec- 
tral resolution, wavelength sampling, and angular resolution, but 
also the simplifications and approximations of the model used to 
. interpret the measurements. 
K^ ■ In this paper we want to evaluate the merits of Milne- 
. ^-H ' Eddington (ME) inversions for the analysis of the polarization 
^^ line profiles emerging from the solar atmosphere. The ME ap- 
?— I I proximat ion does not account for vertical vari ations of the pa- 
. 5^ ! rameters (IUnnoll956l : lRachkovskvlll962L[T967b . so it cannot ac- 
curately describe the solar plasma when rapid changes in height 
are present. What is, then, the significance of the ME parame- 
ters? 

To answer this question it is necessary to simulate the pro- 
cesses of line formation and data inversion. Usually one pre- 
scribes a set of model atmospheres, performs spectral synthe- 
sis calculations, inverts the synthetic profiles, and compares the 
results with the known input. A common approach is to use 
ME models both to generate the spectra an d to invert them 
(e.g.. iNorton et al.ll2006l: iBorrero et"aDl2007h . In that case the 
analysis is internally consistent and the uncertainties of the re- 
trieved ME parameters are mostly due to the noise and, to a 
smaller extent, to the convergence of the algorithm, provided 



that the spectral resolution and wavelength sampling are ap- 
propriate. Uncertainties caused by photon noise are known as 
statistical errors and can be evaluated by means of numeri- 
cal tests or, more efficiently, by using ME response functions 
(lOrozco Suarez & Del Toro Iniestal2007l:lDel Toro Iniesta et al.l 
|201Q|) . However, they represent only a small fraction of the total 
error. Another source of error is the very assumption of height- 
independent parameters, which leads to symmetric line profiles. 
What happens when realistic (i.e., asymmetric) Stokes spectra 
are analyzed in terms of ME models? Do the uncertainties of 
the retrieved parameters increase significantly? Answering these 
questions is the aim of the present work. 

A first study of the cap abilities and limitations of ME in- 
versions was carried out by IWestendorp Plaza et al.l (1 19981) us- 
ing simple (non-ME) model atmospheres. They made a quan- 
titative comparison of results obt ained with the ME code 
of the High Altitude Ob servatory (ISkumanich & LitesI [19871: 
iLites & Skumanichlll99Ql) and t he SIR code (Stokes Inversion 
based on Response functions; iRuiz Cobo & Del Toro Iniestal 
'1992^. The main conclusion of their work was that ME inver- 
sions provide accurate values of the physical parameters aver- 
aged along the line of sight, at least when the stratifications are 
smooth. 



More recently jKhomenko & ColladosI (l2007bl) have investi- 
gated whether the magnetic field stratification itself can be deter- 
mined reliably through inversion of high resolution data. To that 
end, they synthesized the Stokes profiles of the Fe 1 630 nm lines 
with the help of MHD models and inverted them with SIR, al- 
lowing for vertical gradients of the atmospheric parameters. The 
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Fig. 1. Magnetic field strength, inclination, azimuth, and LOS velocity at r = 1 in a simulation snapshot with ifi) - 140 G. Negative 
velocities represent upflows. 



analysis showed that SIR is able to recover the actual magnetic 
stratification for fields as weak as 50 G if no noise is present . 
This work extends the results of IWestendorp Plaza et alJ (Il998l) 
to the case in which the stratifications are not smooth. 

To determine the uncertainties associated with ME inver- 
sions of asymmetric Stokes profiles we use state-of-the-art mag- 
netohydrodynamic simulations (Sect. O. Our goal is to de- 
scribe the solar photosphere as realistically as possible. We 
construct model atmopheres from the simulations and synthe- 
size the Stokes profiles of the Fei 630.2 nm lines emerging 
from them (Sect. O. The SIR code is used for the spectral 
synthesis, so the profiles are asymmetric. Finally, we apply a 
ME inversion to the data (Sect. |4]). In our numerical experi- 
ments, the spatial sampling of the MHD models, 0'.'0287, is 
preserved. There are two reasons why we neglect the eff'ects 
of solar instrumentation: first, they have already been stud- 
ied in the past (e.g., Orozco Suarez et al. 2007, 2010); sec- 
ond, this sampling is close to critical for the observations to 
be delivered by la rge telescopes like t he Advanced Technology 
Solar TelescopeQ (IWagner et al.ll2006l) and the European Solar 



Telescop43 ( |Colladosl [2008|). We invert the pro files with the 
MILOS code (Orozco Sua rez & Del Toro Iniestall2007h Pl A di- 
rect comparison of the retrieved and true parameters allows us 
to determine the eff'ective "heights of formation" of the ME pa- 
rameters (Sect. O and to quantify the errors caused by the ME 
approximation (Sect.©. The conclusions of our work are given 
in Sect. |7l For completeness, the results of ME inversions are 
compared with those of tachogram/magnetogram-like analyses 
in the Appendix. 



2. Magnetohydrodynamic simulations 

We use radiative MHD simulations performed with M URaM, 
the MPS/University of Ch icago RAdiative MHD code (IVogled 
120031: IVogler et alJ l2005h . This code solves the 3D time- 
dependent MHD equations for a compressible and partially ion- 
ized plasma taking into a ccount non-gr ey radiative energy trans- 
port and opacity binning (lNordlund|[l9 82) . 



http://atst.nso.edu/ 



^ http://www.iac.es/project/EST/ 

^ MILOS is programmed in IDL and can be downloaded from our 
website, jhttp://spg.iaa.es/download.asp| 
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Among other problem s, MURaM has b een employed to 
study facular brightenings (JKeller et al.ll2Q04l) . the relation be- 
tween G-band bright points and magnetic flux concentrations 
(ISchtissler et alJ 120031: IShelvag et aP l2004 . the emergence of 
magnetic flux tubes frona the upper convectio n zone to the pho- 
tosphere (ICheung|[20Q6l: ICheung et aDl20Q7l). the strongly in- 
clined magnetic fields of the internetwork ("Sch ussler & Voglej 
I2OO8I). umbral dots ( tSchiissler & Vogl er 200^ solar pores 
(ICameron et alJ l2QQ7l) . and even full sunspots (JRempel et alJ 
l2Q09allbb . MURaM has also been used to evalu ate the diagno s- 
tic potential of spectral lines (Khomenko & Co lladosI I2QQ7 alibi: 
iKhomenko eTaP 12005 alibi: Ishelvag et al.ll2007D . the validity of 
visible lines for th e study of internetwork magn etic fields at high 
spatial resolution (lOrozco Suarez et ^l2007l), and the con tin- 
uum contrast of the solar granulation (iDanilovic et al.ll20Q8l) . 

In this paper we consider a 5 -minute sequence of a mixed- 
polarity simulation run representing a network region with an 
average magnetic field strength (B) = 140 G at log r = -iQxhe 
cadence is 10 s, so we have 30 snapshots. A bipolar distribution 
of vertical fields with (B) = 200 G was used to initialize the 
simulati ons. Additional details ab out this particular run can be 
found in IKhomenko et al.l (l2005al) . 

The computational box has 288 x 288 x 100 grid points 
and covers 6000 km in the horizontal direction and 1400 km 
in the vertical direction. The model extends from z = -800 to 
z = 600 km, with z = km the average of the heights where 
T = 1. The spatial grid sampling is 0^.^0287, implying an equiv- 
alent resolution of 0^.^057 (41.6 km) on the solar surface. The 
simulation provides the density, the linear momentum density 
vector, the total energy density, the magnetic field vector, the 
temperature, and the gas pressure at every grid point. The time- 
averaged radiation flux density leaving the top of the box has the 
solar value F© = 6.34 x 10^^ erg s"^ cm"^. 



3. Spectral synthesis 

In order to compute synthetic Stokes profiles we solve the radia- 
tive transfer equation (RTE) for polarized light. The calculations 
are carried out ii sing the SIR code with the opacity routines of 
IWittmanij (Il974l) . The spectral synthesis is accomplished in two 
steps: first, the input model atmospheres are built from the MHD 
simulations; then, the RTE is solved. 



3.1. MHD models and spectral line synthesis 

The parameters needed for the spectral synthesis are the tem- 
perature, electron pressure, line-of-sight (LOS) velocity, mag- 
netic field strength, inclination and azimuth, and optical depth. 
The simulations provide most of them. However, the electron 
pressure and optical depth need to be computed from the local 
temperature, gas pressure, and density by solving the Saha and 
Boltzmann equations. The optical depth scale is set up assuming 
that z = 600 km (the top of the computational box) corresponds 
to logr = -4.9. This value has been taken fr om the Harvard- 
Smith sonian Reference Atmosphere (HSRA; iGingerich et all 
I1971I) . Finally, the resulting stratifications are interpolated to an 
evenly spaced optical depth grid using second-order polynomi- 
als. The grid extends from logr = -4 to 2 with A logr = 0.05. 
This range of optical depths encompasses the formation region 
of all photospheric lines, except the cores of the strongest ones. 
Figure [T] displays maps of the field strength, inclination, az- 
imuth, and LOS velocity at r = 1 for one simulation snapshot. In 



Li_ 
Q 







Field strength PDF 




10-^ 






- 


10-^ 


: \ 


V 


- 


10-^ 


r 


^^- — ^ 


^: 


in-^ 






, 1 , , 



200 400 600 800 1000 1200 1400 
Magnetic field strength [G] 

Field Inclination PDF 




50 100 150 

Magnetic field inclination [deg] 

Fig. 2. Probability density functions for the magnetic field 
strength and field inclination at log r = -1. 



the velocity map the granulation pattern is clearly visible, with 
granular upflows that are weaker than the intergranular down- 
flows. Some of the small-scale intergranular structures have ve- 
locities of up to 6 km s"^ 

The field strength map shows strong flux concentrations in 
the intergranular lanes. Granules also harbor magnetic fields, but 
they seldom exceed 300 G. There is a tight correlation between 
the field strength and inclination in these simulations: the inter- 
granular fields tend to be vertical, whereas the granules exhibit 
more horizontal fields. Finally, the azimuth map is dominated by 
granular- sized structures with diameters of V'-T' (0.7-1.5 Mm). 

Figure [2] depicts the probability density functions (PDFs J3 
of the magnetic field strength and inclination at optical depth 
logr = -1, averaged over the 30 available snapshots. The field 
strength PDF increases rapidly toward weak fields, indicating 
that most pixels have magnetic fields of the order of hecto- 
gauss. The distribution peaks at about 20 G. The inclination PDF 
shows some vertical fields and a larger occurrence of horizontal 
fields. The simulation run was seeded with mixed-polarity ver- 
tical fields; therefore, the distribution is rather symmetric about 
y = 90°. 



All optical depths refer to the continuum opacity at 500 nm 



^ The PDF is defined such that P(B)dB is the probability of finding a 
magnetic field B in the interval [B, B + dB] . The integral of the PDF is 
unity, i.e., jj" P(5y5 = 1. 
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Table 1. Atomic parameters of the spectral lines. 



Ion 


A (nm) 


Xlow 


logte/) 


Transition 


a 


cr/4 


^eff 


Fei 
Fei 


630.1501 
630.2494 


3.654 
3.686 


-0.75 
-1.236 


5P2 - 5D2 
5Pi - 5Do 


0.243 
0.240 


840.5 
856.8 


1.67 

2.5 



Notes. A stands for the central wavelength of the transition, ;^iow is the excitation potential of the lower atomic level in eV, log gf represents for 
the multiplicity of the lower level times the oscillator strength, a and cr are the collisional broadening parameters in the quantum theory of Anstee, 
Barklem and O'Mara (in units of Bohr's radius, ao), and geff is the effective Lande factor of the transition. 



Once we have constructed model atmospheres for each of 
the 288 X 288 pixels and for all the snapshots, we use them 
to compute the Stokes profiles of Fei 630.15 and 630.25 nm. 
The atomic parameters used in the calculations are given in 
Table [T] Th e lo ggf values have been taken from the VALD 
database (Piskunpv et al. 1995), except for Fe 1 630.25 nm which 
is not available in VALD and comes fro m a fit to the solar spec- 
trum u sing the two-component model of JBorrero & Bellot Rubiol 
(l2002h . The collisional broadening coeflficients a and cr due to 
neutral hydroge n atoms have been e valua ted following t h e pro- 
cedur e of Anstee & O'Maral (Il995h and iBarklem et"^ (11998 ' 
|200Q|) . The abundances have been taken from lTheveninI { 
i.e., a value of 7.46 is employed for iron. 



3.2. Synthesis results 

Figure [3] shows a continuum map of the simulated region. The 
rms contrast, computed as the standard deviation divided by the 
mean, is 14.8% at 630 nm. This contrast exceeds the typical 
values obtained from ground-based observations; the speckle- 
reconstructed G-band images taken at t he Dunn Solar Tele scope, 
for example, show contrasts of 14.1% (lUitenbroek et al.l 12007) . 
Nevertheless, the agreement is reasonable beca use the observed 
values are degraded by instrumental eff^ects (e.g.,lDaniloyic et al.l 



l2008l:IWedemever-Bohm & Rouppe van der \^ortll2009h . 

We have compared the temporally and spatially averaged 
intensity profiles with the corresponding ones in the Fourier 
Transform Spectromete r ( FTS) atlas of the quiet Sun by 
iBrault & NeckeD (Il987h and lNeckell (Il999l) (see Fig.©. For the 
comparison, the simulated profiles have been first normalized to 
unity and then displaced in jvavelength to correct for the solar 
gravitational redshift (see Lppresto et al. 1980). Also, an addi- 
tional minor correction to the wavelength shift has been allowed 
to improve the fits. The figure shows that the widths of the syn- 
thetic profiles are very similar to those measured in the FTS at- 
las. However, the bottom panels indicate that the simulations do 
not completely reproduce the line asymmetries: the diff'erences 
between observed and synthetic profiles are small (of the order 
of 2.5%), but not zero. 

In summary, despite some diff'erences between the FTS and 
the synthetic profiles, the simulations appear to explain the ob- 
servations rather well. Therefore, we hope that they also pro- 
duce sufficiently realistic polarization spectra, in particular with 
regard to their asymmetries. 
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Fig. 3. Continuum intensity map for the simulation snapshot de- 
picted in Fig.[T] The wavelength is 630 nm. 
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Fig. 4. Average intensity profiles from the simulations (solid) 
and the FTS atlas (dashed), for the Fei lines at 630.2 nm. The 
bottom panel shows the intensity diflferences (FTS - simulation) 
in percent. 



4. ME inversion of the Stokes profiles 

To determine the vector magnetic field and the LOS velocity 
we perform ME inversions. Given the high spatial and tempo- 
ral resolution of the simulations, the macroturbulence is set to 
zero and the filling factor to unity, i.e., the magnetic atmosphere 
is assumed to occupy the whole pixel (one-component model 
atmospheres). A total of 9 quantities are determined from the in- 



version: the thermodynamic parameters So, Si, j]o, A Ad, and a 
(representing the intercept and the slope of the source function, 
the line-to-continuum opacity ratio, the Doppler width, and the 
damping parameter), the strength, inclination, and azimuth of the 
magnetic field vector (B, y, and;^), and the line-of-sight velocity 
(vlos)- The Stokes profiles are taken from a single snapshot of 
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the simulation. No noise is added to the spectra to better isolate 
the uncertainties due to the ME approximation. 

The Fe i lines at 630 nm belong to the same multiplet and are 
formed under very similar thermodynamic conditions. Hence, 
we can reliably assume that the ME thermodynamic parameters 
are the same for both lines except for 7/0 . A simultaneous in- 
version with no extra free parameters is thus possible using a 
constant rjo ratio and the same AAd and a values for the two 
lines. This is the strategy implemented i n several ME codes (for 
details, see lOrozco Suarez et al. 11201 Obi) . 

We use the same initial guess model for all the pixels, stop- 
ping the inversion when convergence is achieved or 200 it- 
erations have been performed. The initial model is given by 
^0 = 0.2, ^1 = 0.8, 7/0 = 6.5, A.1d = 30 mA, a = 0.03, 
5 = 200 G, r = 20°, ;t = 20°, and vlos = 0.25 kms'^ 

5. Understanding ME inferences 

The line asymmetries observed in the solar atmosphere imply 
that there are vertical gradients of the physical parameters, be- 
cause the line wings trace relatively deep layers and the core is 
formed higher in the atmosphere. By contrast, the ME fits are 
strictly symmetric and deliver height-independent atmospheric 
parameters. Thus, it is important to verify that the ME model 
can indeed be used for the analysis of observations at very high 
spatial resolution. 

Figure \5\ shows the magnetic field strength, inclination, az- 
imuth, and LOS velocity stratifications in three pixels of the 
simulation snapshot, labeled (a), (b), and (c). The figure also dis- 
plays the corresponding Stokes /, Q, U, and V profiles. The re- 
sults of the ME inversions are overplotted in red. Case (a) shows 
symmetric Stokes profiles, in (b) the profiles are rather asymmet- 
ric, and (c) shows three-lobed Stokes V spectra together with 
anomalous linear polarization signals, (a) represents a strong 
field case and (b) and (c) correspond to weak fields. In the three 
examples the atmospheric quantities undergo large variations 
with optical depth. 

The ME fit is good in (a) and worse in (b) and (c). Clearly, as 
the asymmetry level increases, the ME model has more difficul- 
ties to reproduce the profiles. The misfits are obvious in Stokes 
2, U, and y, and less dramatic in Stokes /. 

The models retrieved from the inversion are shown in the 
left panels of Fig. [5] (red lines). The height-independent ME 
parameters can be interpreted as averages of the actual strat- 
ifications weighted by the co rresponding response functions 
(IWestendorp Plaza et al.lll998l) . In general it is difficult to con- 
firm this by simply looking at the atmospheric stratifications, but 
case (c) provides a particularly clear example. This case repre- 
sents a pixel whose ME fit is not satisfactory. The analysis of the 
stratifications shows that the profiles arise from an atmosphere 
that has sharp discontinuities in field strength, inclination, az- 
imuth, and LOS velocity, located more or less at the same op- 
tical depth. The two parts of the atmosphere separated by the 
discontinuity leave clear signatures in the emergent Stokes V 
spectra, to the point that the magnetic and kinematic properties 
of the plasma above and below the discontinuity can roughly be 
guessed from a simple inspection of the profiles: there is a weak 
field associated with small velocites and a stronger field showing 
large redshifts. Surprisingly, however, the ME model returned by 
the inversion seems to describe only the weak field. 

To understand why this happens, it is important to realize that 
the inversion algorithm uses all the wavelength samples to deter- 
mine the best-fit ME parameters. As mentioned before, diff'erent 
wavelength positions across the line trace diff'erent atmospheric 



layers. Thus, the ME inversion is forced to return average pa- 
rameters along the LOS in order to fit the whole line profile rea- 
sonably well without any bias toward better fits in the core or the 
wings. This favors the weak field component of the atmosphere 
because it occupies most of the line-forming region in this par- 
ticular example. 

Figure [5] demonstrates that the ME model parameters coin- 
cide with the real stratifications at specific optical depths. This 
allows us to define the eff'ective "height of formation" of the ME 
parameters. 

Formation-height maps have been calculated by taking the 
optical depth at which the stratification is closer to the inferred 
ME parameter. The formation heights are constrained to be in the 
range from log r = to -2 because this interval includes most of 
the layers to which the Fe i lines are sensitive. When more than 
one value of the MHD stratification coincides with the corre- 
sponding ME parameter, we select the one located deeper in the 
atmosphere. The optical depth of the minimum (or maximum) 
of the MHD stratification is taken if the ME parameter is smaller 
(or larger) than all stratification values. Even in that case, the for- 
mation height is not allowed to go outside of the interval from 
logr = to -2. 

Figure[6] shows the results for the magnetic field strength and 
the LOS velocity. For convenience, we also display the contin- 
uum image and the field strengths retrieved from the inversion. 
Diff'erent colors indicate diff'erent atmospheric layers. There are 
clear diff'erences between the two maps: the granular centers are 
predominantly green (logr e [-0.7, -1.2]) in the velocity map 
and green-red (logr g [-0.3, -0.7]) in the field strength map, 
demonstrating that the ME inversion extracts the velocities from 
higher optical depths than the magnetic fields, at least in gran- 
ules. The intergranular lanes tend to show blue colors in the two 
maps (logr e [-L3,-L7]). The sharp red-yellow discontinu- 
ities observed at the border of granules in the velocity map have 
little meaning since they are a result of the way we estimate for- 
mation heights. These regions have zero velocities (Fig. [T]) and 
also zero velocity gradients along the LOS, which makes our al- 
gorithm to select layers close to the bottom of the photosphere 
(limited to logr = 0). Both maps exhibit significant pixel-to- 
pixel diff'erences, especially the field strength map. This "noise" 
is due to MHD stratifications with many jumps in the vertical 
direction. 

In summary, ME inversions provide results that can- 
not be assigned to a constant optical depth. The same pa- 
rameter may show formation-height diff'erences of up to 
1-1.5 dex across the FOV. Also, the heights to which 
the ME parameters refer change depending on the par ame- 
ter, as predicted by iDel Toro Iniesta & Ruiz Cobol (Il996h and 
ISanchez Almeida et al.l (Il996h . In the case of the Fei 630.2 nm 
lines, we find mean optical depths of logr = -1.0 and -1.1 for 
the LOS velocity and the field strength, respectively. This in- 
cludes granular and intergranular regions. If only intergranular 
regions are considered, the mean optical depths move -0.2 dex 
toward higher layers. The rms variation of the formation heights 
is 0.4 and 0.5 dex, respectively. 



6. Inversion results 

In what follows we compare the ME inversion results with the 
original MHD models. To that end we use the atmospheric pa- 
rameters of the simulations at log r = -I. This layer corresponds 
to the average formation height of the ME parameters. As such, 
it represents the best choice for a "reference model". 
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Fig. 5. Examples of MHD atmospheres and simulated profiles (black) and ME fits (red) for three diff'erent pixels. Left panels: 
magnetic field strength, inclination, azimuth, and LOS velocity stratifications. The red horizontal lines indicate the inversion results. 
Right panels: Stokes /, 2, U and V profiles synthesized from the MHD simulations with no noise (black) and ME fits (red). Cases 
(a), (b) and (c) correspond to (x,y)=(L75, LOS), (2.90, 2.58), and (L67, L73) Mm in Figs. [6] and El 



Figure [7] displays maps of the magnetic field strength, incli- 
nation, azimuth, and LOS velocity in the reference model (left 
column) and the models retrieved from the inversion (right col- 
umn). To better visualize the details we only show a small area 
of about 9 Mm^. The strong resemblance between the reference 
parameters and the ME models is obvious: the shapes of the dif- 
ferent structures are well reproduced and only small diff'erences 
can be recognized. Sometimes the inversion yields bad results 
for the inclination and azimuth, but this happens mainly in areas 
with weak polarization signals. 

From a visual inspection of the maps, one can say that the 
ME inversion is able to determine the magnetic field vector satis- 



factorily. Even structures with field strengths as low as 100 G are 
well recovered. To make more precise statements. Fig. [8] shows 
the parameters inferred from the fit vs the MHD parameters at 
log T = - L These scatter plots allow us to estimate the uncer- 
tainties that can be expected from the use of the ME approxima- 
tion, since no noise has been added to the profiles (Sect.©. 

As can be seen, the scatter is larger for the magnetic field in- 
clination and azimuth. The mean values of the parameter^ (red 
dots) show that the magnetic field strength is really close to that 
in the reference model from to 500 G. For stronger fields the 



^ The average values have been calculated by taking bins along the X- 
axis of size 28 G, 3°, and 115 m s~^ , depending on the physical quantity. 
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Fig. 6. Top: Optical depths at which the inferred ME parameters coincide with the real magnetic field strength and LOS veloc- 
ity stratifications (left and right, respectively). Bottom: Magnetic field strengths retrieved from the ME inversion and normalized 
continuum intensities (left and right, respectively). 



retrieved values are slightly underestimated, although the devia- 
tion does not exceed ~250 G. The red lines indicate maximum 
rms fluctuations of about ~80 G in the whole range of strengths. 
The inclination shows rms fluctuations smaller than 10° for ver- 
tical fields and 15° for inclined fields. The rms variation in the 
azimuth is about 15°. The LOS velocity panel shows that the 
retrieved velocity is some 200-300 m s"^ smaller than the refer- 
ence velocities for receding flows (intergranular lanes). The rms 
values are smaller than ~ 500 m s"Mn the full velocity range. 

The scatter observed in the various panels of Fig.[8]originates 
from the use of ME model atmospheres (unable to fit asymmetric 
Stokes profiles) and the pixel-to-pixel variations of the forma- 
tion height of the ME parameters, as explained in the previous 
section. The deviation of the ME field strengths from a one-to- 
one correspondence with the MHD models can easily be under- 
stood by looking at the top panel of Fig. [6] and recalling that we 
have chosen the atmospheric layer at logr = -1 as a reference. 
In those spatial locations where the optical depth assigned to 
the retrieved ME parameter is smaller than the optical depth of 
the reference layer, the resulting field strength will "apparently" 
be underestimated. These spatial locations are associated with 
strong flux concentrations; in the MHD models they spread out 
with height, therefore we retrieve weaker fields. 



The rms diff'erences between the inversion results and the 
MHD models tell us how much a single ME parameter could 
deviate from the real value, but only if the mean differences are 
zero or close to zero. Non-zero mean difl'erences indicate that 
systematic errors exist. For this reson, when the mean difl'erence 
is larger than the standard deviation, the former should be pre- 
ferred as a better estimate of the true error. 

The choice of log r = - 1 for the reference model is appro- 
priate because it produces the smaller mean and rms values on 
average. To illustrate this, Fig.[9]represents histograms of the dif- 
ferences between the inferred parameters and the MHD model at 
three optical depths (logr = -0.5, -1,-1.5, coded in black, red, 
and green, respectively). 

For the magnetic field strength, the histogram corresponding 
to logr = -1 peaks around zero. The maximum shifts toward 
negative values when the inversion results are compared with 
deeper layers (fields are underestimated on average) and toward 
positive values when the comparison is made with higher lay- 
ers (over-estimating the strength). The full width at half maxi- 
mum (FWHM) of the distribution is about 30 G for log r = -1, 
and increases up to -45 and -50 G for logr = -1.5 and -0.5, 
respectively. These eff'ects are less pronounced for the field in- 
clination: the peaks of the histograms are located at zero and 
the FWHM varies from 6° (logr = -1.5) to -13° and -23° 
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Fig. 7. From top to bottom: magnetic field strength, inclination, 
azimuth, and LOS velocity. The left column represents the MHD 
parameters at log r = -I. The right column shows the results of 
the ME inversion of the Fe i lines at 630.2 nm. 



(logr = -1 and -0.5, respectively). The larger FWHMs origi- 
nate from the extended wings of the distributions. The azimuth 
histogram does not vary when the comparison is made with dif- 
ferent optical depths. In this case, the FWHM is about 20°. 

The histograms of the LOS velocity diff'erences show larger 
variations. The one corresponding to logr = -1 has the smaller 
FWHM (~ 500 m s"^). It also features a long tail toward negative 
values caused by pixels located in intergranular lanes. The asym- 
metry of the histograms around the location of the peaks changes 
dramatically when we compare the inversion results with dif- 
ferent atmospheric layers. For instance, if the reference layer 
is taken at logr = -0.5, the histogram is a clear combination 
of two diff'erent distributions, one representing granular centers 
(higher and narrower) and the other representing intergranular 
lanes (smaller in amplitude and broader). 



7. Summary and conclusions 

In this paper we have analyzed radiative MHD simulations of the 
quiet Sun. We have used them to synthesize Stokes line profiles 
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Fig. 8. Scatter plots of the magnetic field strength, inclination, 
azimuth and LOS velocity inferred from the ME inversion vs 
the MHD parameters at log r = -I. The dashed lines represent 
one-to-one correspondences. The red dots are the mean values 
of the parameters in small, evenly-spaced bins along the X-axis. 
The blue lines represent the rms fluctuations of the parameters 
in those bins. 
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Fig. 9. Normalized histograms of the difl'erences between the in- 
ferred ME model parameters and the real ones taken at difl'erent 
optical depths. 



in three difl'erent spectral regions (525.0, 617.3, and 630.2 nm). 
The comparison of the synthetic profiles with the FTS atlas sug- 
gests that the simulations describe quite satisfactorily the physi- 
cal conditions of the solar photosphere, although the MHD mod- 
els are slightly hotter than the HSRA around r = L 

After synthesizing the Stokes profiles, the applicability of 
ME inversions to high spatial resolution observations has been 
examined. We have considered the case of the Fei lines at 
630.2 nm. The analysis of the profiles by means of ME inver- 
sions has allowed us to characterize the uncertainties that can be 
expected from the ME approximation. For this reason, the syn- 
thetic profiles were not degraded by noise, instrumental eff'ects, 
or spatial resolution. 

The main limitation of ME inversions is that they provide 
constant atmospheric parameters, whereas the MHD models fea- 
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ture physical properties that change with height. This limitation 
means that ME models are unable to reproduce spectral line 
asymmetries. Consequently, the ME inferences cannot be as- 
signed to a specific optical depth. Depending on the conditions 
of the atmosphere, the retrieved ME parameters sample diff'erent 
layers. 

However, from a statistical point of view we conclude that 
ME inversions provide fair estimates of the physical conditions 
prevailing at log r ~ - 1 . The rms uncertainty is smaller than 
30 G for the magnetic field strength, 13° and 20° for the field in- 
clination and azimuth, and 500 m s"^ for the LOS velocity. Thus, 
ME inversions are appropriate for statistical analyses of the so- 
lar photosphere. This being said, it is important to realize that the 
errors may be large for individual pixels, even if the best-fit pro- 
files reproduce the observations satisfactorily (the field strength 
in case a of Fig. [5] is a good example of this). 

Finally, we want to stress that the uncertainties associated 
with t he ME approximation are larger than those due to photo n 
noise (lOrozco Suarez et al. ll2Q06l: iDel Toro Iniesta et aLlboiOl) . 
However, the noise has another undesirable eff'ect: it hides the 
weaker polarization signals. This fact has not been considered in 
our study. 

Appendix A: Milne-Eddington vs classical proxies 

The study of the solar atmosphere relies on the availability of 
precise magnetic fields and LOS velocities. Thus, one needs ro- 
bust diagnostics in order to extract this information from the 
Stokes spectra. Classical methods such as tachogram techniques, 
the weak field approximation, and the center-of-gravity tech- 
nique represent an alternative to Stokes inversions. 

For some of these methods, the random errors induced by 
photon noise have been estimated not to exceed ~ 20 m s"Mn 
the case of the LOS velocity or ~ 10 G in the case of the mag- 
netic flux (see e.g., Scherrer et al. 1995; Scherrer & SDO/HMI 
Team 2002; Martinez Pillet 2007). However, like in the case 
of ME inversions, systematic uncertainties coming from the hy- 
potheses underlying the technique are expected to be larger than 
the random errors themselves. A thorough study, similar to the 
one we have performed for the ME technique, is thus in order. 
We carry out such an analysis in this Appendix for the Fourier 
tachometer technique (Beckers & Brown 1978; .Brownlfl9 8lLl, 
the center-of-gravity method ( Semelll967l:lRees & Sem el 1979), 
and the weak field approximation ([Landi degl' Innocent! .1992: 
iLandi Degrinnocenti & Landolfill2004l) . 

Both the center-of-gravity method and the weak field ap- 
proximation are applied to the whole profiles while the Fourier 
tachometer uses only four wavelength samples across the Stokes 
/ profile (-9, -3, 3, and 9 pm). The center-of-gravity technique 
extracts the longitudinal component of the magnetic field, ^z, 
from the separation between the barycenters of the Stokes I -\-V 
and I - V profiles, ^z can also be obtained with the weak field 
approximation through a proportionality between the Stokes V 
profile and the wavelength derivative of Stokes /. The transverse 
component of the field, in turn, is derived through a proportion- 
ality between Stokes L and the second wavelength derivative of 
Stokes /. Regression fits are used between the circular (lin- 
ear) polarization profiles and the first (second) derivatives of the 
intensity profiles for increased accuracy. Then, the magnetic in- 
clination is obtained from the ratio between the transverse and 
longitudinal components of the field. 



^ We in fact use the formula proposed bv lFerna ndes' (*1992|). 
^ Stokes L is the total linear polarization, L - -^Q^ + V^. 



Figure [TOl summarizes the results. Each column refer to a 
diff'erent set of lines: Fe i 630.15 nm (left), Fe i 630.25 nm (mid- 
dle), and the two lines simultaneously considered (i.e., ME in- 
version; right). The labels on the ordinates are self-explanatory, 
while the abscissae give the values of the corresponding quanti- 
ties at log T = - 1 . 

The less accurate method turns our to be the weak field ap- 
proximation: the inferred magnetic inclinations show larger rms 
fluctuations than the ME ones, and the longitudinal component 
of the field displays a clear saturation for fields stronger than 
1000 - 1100 G when calculated with the line at 630.25 nm. 
For ^z values above 1 kG, the weak field inferences resulting 
from the 630.15 nm line seem to be closer to the MHD values at 
log T = -1 than the ME ones. Nevertheless, they present a larger 
scatter. The center-of-gravity method looks very robust and, in- 
deed, it presents less scatter than ME inversions in the case of 
the longitudinal field component and the LOS velocity (only the 
results from the 630.15 nm line are shown; the results for 630.15 
nm are very similar). The good perfo rmance of th e cente r-of- 
gravity method was noticed earlier by ICauzzi et"aD (1 19931) and 
Uitenbroek (2003). Unfortunately, this technique does not pro- 
vide information about the field inclination. The LOS velocities 
resulting from the tachometer are fairly comparable to those of 
the ME inversion. The scatter is similar in both cases. 

In summary, the ME inversion seems to be the more com- 
plete and accurate technique although none of the others can be 
discarded. In particular, a combination of the center-of-gravity 
technique for calculating B^ and vlos along with the weak field 
approximation for the magnetic inclination may represent a suit- 
able alternative which is much less expensive in terms of com- 
puting resources. It is important to note, however, that the results 
of our study are only valid when the magnetic field is spatially 
resolved. Further investigation is needed to check the applicabil- 
ity of these techniques when the field is unresolved. This addi- 
tional investigatio n is important in view of the theoretical d evia- 
tions predicted by lLandi Degl'Innocenti & Landolfil (|2004|) . 

Acknowledgements. This work has been funded by the Spanish MICINN 
through projects AYA2009-14105-C06-06 (including European FEDER funds) 
and PCI2006-A7-0624, by Junta de Andalucia through project P07-TEP-2687, 
and by the Japan Society for the Promotion of Science. 



References 

Anstee, S. D., & O'Mara, B. J. 1995, MNRAS, 276, 859 

Barklem, R S., Anstee, S. D., & O'Mara, B. J. 1998, PASA, 15, 336 

Barklem, R S., Piskunov, N., & O'Mara, B. J. 2000, A&AS, 142, 467 

Beckers, J. M., & Brown, T. M. 1978, Future solar optical observations needs 
and constraints, 189 

Borrero, J. M., & Bellot Rubio, L. R. 2002, A&A, 385, 1056 

Borrero, J. M., Tomczyk, S., Norton, A., Darnell, T., Schou, J., Scherrer, P., Bush, 
R., & Liu, Y. 2007, Sol. Phys., 240, 177 

Brault J., Neckel H., 1987, Spectral atlas of solar absolute disk- 
averaged and disk-center intensity from 3290 to 125 lOA, ftp.hs.uni- 
hamburg.de/pub/outgoing/FTS-atlas, (Hamburg: Hamburg Univ.) 

Brown, T. 1981, in Solar instrumentation: What's next?, ed. R.B. Dunn, (NSO: 
Sunspot, NM), 150 

Cameron, R., Schussler, M., Vogler, A., & Zakharov, V. 2007, A&A, 474, 261 

Cauzzi, G., Smaldone, L. A., Balasubramaniam, K. S., & Keil, S. L. 1993, 
Sol. Phys., 146, 207 

Cheung, M. C. M. 2006, PhD Thesis, University of Gottingen, Germany, 
http ://webdoc . sub. g wdg . de/dis s/2006/cheung/ 

Cheung, M. C. M., Schussler, M., & Moreno-Insertis, F. 2007, A&A, 467, 703 

CoUados, M. 2008, Proc. SPIE, 7012, 70120J 

Danilovic, S., Gandorfer, A., Lagg, A., Schussler, M., Solanki, S. K., Vogler, A., 
Katsukawa, Y., & Tsuneta, S. 2008, A&A, 484, L17 

Del Toro Iniesta, J. C. 2003, Introduction to Spectropolarimetry. Cambridge, 
UK: Cambridge University Press 

Del Toro Iniesta, J. C, & Ruiz Cobo, B. 1996, Sol. Phys., 164, 169 



Orozco Suarez et al.: Milne-Eddington inversions of high resolution observations of the quiet Sun 



630 J 5 nm 



630.25 nm 









Both 


lines 




„ 2 




a = 80 G 




y 


C 1 








^f^ 


o 1 








4g^ 


'^ 






^ 


W^ 


UJ 

2 






1^ 




|-. 




-"'? 






^ -2. 


^ 


y 







-1 1 

Bz ot log T = -1 [kG] 



-2 



70 G 



a = 70 G 




-1 1 

Bz ot log T = -1 [kG] 



o 
o 
u 

E 
o 



o 

5 








3 




-2 -1 
Velocity 


1 
ot log T = 


2 3 
-1 [km/s] 


4 


> 


4 
3 
2 

1 


-1 

-2 

-3 




a 


= 330 


m/s 


^ 


^= 


E 


■ i^rfn^ 


^^; 


S 




>N 


M 


PP^ 


1 


^ 





~ 


3 


-2 -1 
Velocity ot log 


1 

T = 


2 3 4 
-1 [km/s] 


150 


1 


a = 14" 




^ 






■ li^K^ 






. 


.driiHty^ 




"*' ^' 


jdi 


iP^- 


100 


f-- -'1^ 


P 


y^F 


50 


IF^ 


^ ^ 




/i^^^ 




!^.- 


0, 


1^ -- 







50 100 

Inclinotlon ot log t 



50 100 

Inclinotlon ot log t 



150 
CI 



Fig. 10. Scatter plots of the longitudinal magnetic field, B^ (first two rows), of the LOS velocity (third row), and of the magnetic 
field inclination (fourth row), as functions of their corresponding values in the MHD simulations at logr = -1. The first column 
refers to results obtained with the 630.15 nm line; the second column displays results obtained with the 630.25 nm line; results with 
the two lines at the same time are shown in the third column. Color and line codes are the same as in Fig.O The mean value of the 
rms diff'erence between the inferences and the simulations at log r = - 1 is given in the upper left corner of each panel. 
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